Vortex- Antivortex Lattice in Ultra-Cold Fermi Gases 
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We discuss ultra-cold Fermi gases in two dimensions, which could be realized in a strongly con- 
fining one-dimensional optical lattice. We obtain the temperature versus effective interaction phase 
diagram for an s-wave superfluid and show that, below a certain critical temperature T c , spontaneous 
vortex-antivortex pairs appear for all coupling strengths. In addition, we show that the evolution 
from weak to strong coupling is smooth, and that the system forms a square vortex-antivortex lattice 
at a lower critical temperature Tm- 

PACS numbers: 03.75.Ss, 05.30.Fk, 34.50.-s, 32.80. Pj 
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The presence of quantized vortices is a strong indi- 
cation of the existence of a superfluid state. The re- 
cent experiment from the MIT group marked the first 
observation of vortices in atomic Fermi gases. This re- 
sult complements previous experiments involving s-wave 
ultra-cold Fermi gases 0, H0- These studies 

combined Q, 0, 0> M, M> correspond to the exper- 
imental realization of the theoretically proposed BCS- 
to-BEC (weak-to-strong coupling) crossover in three di- 
mensional continuum s-wave superfluids HS0. More 
recent generalizations of these earlier theoretical re- 
sults indicate that the critical temperature in a trap ge- 
ometry is higher than in the continuum case through- 
out the BCS-to-BEC evolution. S-wave Fermi systems 
have also been investigated in one dimensional (ID) [l2^ 
and three-dimensional (3D) 13] optical lattices. Further- 
more, the effects of dimensionality have been analyzed in 
p-wave Fermi and Bose |l5| systems. In the par- 
ticular case of Bose systems, the presence of topological 
defects associated with the phase of the order parameter 
has been detected in a nearly two-dimensional geometry 
[l5j . These optical lattice experiments 0,0,0,0 are 
now allowing the exploration of dimensional (ID, 2D, and 
3D) and correlation effects in interacting Fermi gases. 

Given all the recent advances in experimental tech- 
niques, we discuss the evolution from weak to strong cou- 
pling superfluidity of confined 2D s-wave ultra-cold Fermi 
gases. We show that a Berezinksii-Kosterlitz-Thouless 
(BKT) transition occurs at finite temperatures, 
and that the strong coupling limit produces a critical 
temperature Tbkt = 0.5ei?, where cf is the Fermi en- 
ergy. In the strong coupling limit of a 2D Fermi gas, 
the superfluid transition is not characterized by Bose- 
Einstein condensation (BEC) as in 3D, but by the BKT 
transition. Below Tbkt, pairs of vortices and antivortices 
appear spontaneously for all couplings, and eventually 
condense into a vortex-antivortex (VA) square lattice as 
the temperature is lowered further. The lattice melting 
temperature is shown to be Tm ~ 0.1 cf in the strong 
coupling limit, and the melting mechanism is controlled 
by dislocations 000. 

The appearance of spontaneously generated VA bound 



states and the existence of the VA square lattice are ma- 
jor characteristics of the 2D physics proposed here. This 
spectacular effect could be measured via density or veloc- 
ity sensitive techniques. Vortices or antivortices could be 
detected in a density sensitive experiment without stir- 
ring the condensate, but the topological charge associ- 
ated with the sense of rotation should be detectable only 
in velocity sensitive experiments. 

Two-dimensional Fermi systems can be prepared by 
means of a ID optical lattice, where tunnelling between 
lattice sites is suppressed by a large trapping potential. 
A typical trapping potential is Vtrap = —Vq exp[— 2(x 2 + 
y 2 )/w 2 ]cos 2 (k z z), where A z = 2ir/k z is the wavelength 
of the light used in the laser beam. We assume that the 
width w is such that w ^ Xf, where \p = 2n jkp is pro- 
portional to the interparticle spacing of a Fermi gas with 
Fermi wavevector k^?. The atom transfer energy along 
the z direction can be estimated using a simple WKB 



expression leading to t z ss Wo exp 



W )/E r 



where W = fuv /2 = ^JV Q E r /2. Here, E r = h 2 k 2 /2M 
is the recoil energy of atoms with mass M. Thus, in the 
2D limit of t z — > 0, there is no finite critical temperature 
for BEC, and the superfluid phase in both the BEC and 
BCS regimes is characterized by VA bound pairs. 

To explore the physics described above, we consider 
a 2D continuum model of Fermi atoms of mass M and 
density n = kp/2ir, with Hamiltonian (H = ks = 1) 



Vkk'&k q 6k'q, 



(1) 
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where & kq = V ! -k+q/2,4.V'k-i-q/2,t and £ k = «k - M) witn 
energy dispersion e k = k 2 /2m and chemical potential ul. 
Following the procedure discussed in detail in Ref. |20| . 
we obtain the following separable expression for the in- 
terparticle potential in fc-space, 



Vkk< = -Ar(k)r(k'), 



(2) 



where A is the interaction strength and, in the particular 
case of s-wave symmetry (i = 0), r(k) = (1 + fc/fco) , 
with i?o ~ k^ 1 pl a yi n g the role of the interaction range. 
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The partition function Z at temperature T = f3 1 
is an imaginary-time functional integral with action 
S = /jf dT^2 Ka ^ a {T)d T ^ a {T) +H}. Introducing the 
Hubbard-Stratonovich field q (r), which couples to 
ip'ip', and integrating out the fermionic degrees of free- 
dom, we obtain 



Z = J V<j>V<F exp(-S eff [<M*]), 



with the effective action given by 
rf> — '\Mt)\ 2 

k 



f 

Sc.s = dT /~2 
Jo u 



A 



& -Tr mG k Ur) 



The symbol Tr denotes the trace over momentum, imag- 
inary time and Nambu indices, and G^ k ,(r) is the (in- 
verse) Nambu matrix, 



-(d T +Ck)(5k,k' A k ^ k /(r) 



-(9 T - £k)<5k,k' 



(4) 



with A k , k ,(r) = k _ k ,(r)r((k + k')/2). 

Assuming q (r) = Ao<5 qj o + ?7 q (T) and performing an 
expansion in S e s to quadratic order in 77, one can write 
the effective action as S e g = Sq[Aq) + Sa uc t[r], rj*], where 
5o[A ] is the saddle point action, and the fluctuation 
action Sfluctfo, V*} can be expressed in terms of ampli- 
tudes and phases via r](q) = \r](q)\e lS ( q \ Integrating out 
the amplitudes, and Fourier transforming to position and 
imaginary time r = (r, r), one obtains 

Sfiuct = \ f dr [ Pij did^djOir) + A{d T df - iBd T 9(r)] . 



In this expression, 
A(n,A ,T) = -^J2 



^tanhf^ 
El \ 2T 



El 



where Y k = (2T)- 1 sech 2 (E w /2T) is the Yoshida distri- 
bution, and 

B ( ( ,.A ,T) = ±£| to »h(§ 



The quantity p^ represents the superfluid density tensor 
Pij (fi, A , T) = -L £ [2n Q (k)d i d j Z k - Y^di^d^] , 



where di denotes the partial derivative with respect to 



hi, and 



n (k) 



1 



tanh ( — ^ 
E k \2T 



is the momentum distribution. Notice that p xx = Pyy 
p s , while p xy = p yx = 0. 



The decomposition of 6{t,t) = v (r) + 9 s (t,t) into 
a static vortex part 9 v (r) and a spin- wave part O s (t,t) 
permits us to rewrite 



where S„ 



Sauct — S v + S 

s w 

\jdrp s [W6 v (r)Y 



(5) 



2 

The 



while S sw = 
(3) hi dr [Ps [V^(r)] 2 + A [d T 8 sw (r)] 2 - iBd T 9 sw {r) . 

spin-wave part can be integrated out to give 
= Eqrin[l-exp(-u; q /T)], where w(q) = c|q| 
is the frequency and c = <J p s j A is the speed of the 
spin-wave. Here, Q sw is the spin-wave contribution to 
the thermodynamic potential f2 = f2 + £l sw + &v i where 
fio = TSq[Aq] and fl v are the saddle-point and vortex 
parts, respectively. 



The self-consistent equations for p, Ao and T c can be 
derived from the effective action S c g = So + S sw + S v 
as follows. The order parameter equation is obtained 
through the stationarity condition [SS e s/S<p^(T , )]A = 0, 
leading to 



A ^ 



HWtanhf^ 
2£ k \ 2T 



(6) 



where i5 k = \/£ k + |A k | 2 is the quasiparticle excitation 
energy, and A k = Aor(k) plays the role of the order pa- 
rameter function. Elimination of the interaction strength 
A in favor of the two-body bound state energy E^ in vac- 
uum is possible through the relation 



i = v |r(k)| 2 

A 4" 2e k - E b ' 

k 



(7) 



The number equation is obtained via N p = —dfl p /dp, 
where ft p = CIq + fl sw , leading to 

N p = N + N sw . (8) 

Here, N = -dQo/dfi = 2Y^ k no(k), and N sw = 
-dVt sw /dp. 

The equation for the critical temperature T c = Tbkt 
is determined by the Kosterlitz-Thouless |17| condition 



Tbkt — -^Ps(p, Aq, Tbkt)- 



(9) 



The self-consistent solutions of Eqs. iJBJ , JSJ and © de- 
termine /i, Ao and Tbkt as functions of the two-body 
binding energy E\, (or interaction strength A). Solu- 
tions for Tbkt are shown in Fig. ^ The curves labeled 
MF represent the mean-field (saddle-point) solution for 
T c = Tmf, which is obtained by solving only Eqs. © and 
© with A = 0. The disparity between the BKT and 
MF solutions is larger with increasing A (or |-E&|), indicat- 
ing that strong coupling results are dramatically affected 
by phase fluctuations. Furthermore, Tbkt = 0.5ef in the 
strong coupling limit where p ~ n/M, with n = Nq/L 2 
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FIG. 1: Mean field (MF), Berezinskii-Kosterlitz-Thouless 
(BKT) and vortex-antivortex lattice melting (M) critical tem- 
peratures (in units of ep) as functions of (a) absolute value 
of the binding energy \Eb\ (in units of ep) and (b) interaction 
strength A (in units of g^, where <?2D is the two-dimensional 
density of states). 
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FIG. 2: (a) Chemical potential and (b) order parameter am- 
plitude (both in units of ep) evaluated at the transition tem- 
perature Tbkt, as functions of \Eb\ (main figure) and of the 
interaction strength A (in units of g^u) (inset). 



being the fermion density. This shows that although 
there is no finite T c for BEC in 2D, there is still a su- 
pcrfhiid transition with a high T c oc n. Finally, ^(Tbkt) 
and Ao(Tbkt) are shown in Fig.|21as functions of \Ef,\ or 
A. Notice that ^(Tbkt) (Ao(Tbkt)) is a monotonically 
decreasing (increasing) function of \Eb\ or A, and that fi 
changes sign at \E b \ » 2.23 (A » 0.419). 

Since the vortex part is directly related to the trans- 
verse velocity Vt(r) = V#„(r), where V • ft(r) = 0, we 
can express the vortex contribution to the action via 
V x v t (r) = 27rz J2i TiiS(r — ri), where n, = ±1 is the vor- 
tex topological charge (vorticity). Thus, we focus next 
only on the vortex action, which becomes 

Sv = y~ = (2^)f E WjGfr *•;■) + E E * n l 

where H v is the vortex Hamiltonian, E c is the vortex core 
energy, and G(rj — r'-) plays the role of the interaction 
potential between topological charges rii and rij , and sat- 
isfies VjG(r-r') = 6(r — r'). Given that in the vicinity of 
Tbkt VA pairs are formed, it is possible that they crystal- 
ize into a solid at a lower temperature Tm- We note that 
the vortex action just obtained allows for the appearance 



of a VA solid (with well defined lattice structure) at a sta- 
ble minimum of the interaction potential. At T <C Tm, 
the lattice configuration that is compatible with short 
range core repulsions is such that vortices and antivor- 
tices form two square sublattices of side b = ay2, where 
2o is the VA pair size. This VA lattice can be constructed 
as a superposition 9l(x, y) = JV j CijOvAix—ib, y—jb/2), 
where 

2ay 



9y\(x,y) = arctan 



IT 



For convenience, we choose the line connecting the vortex 
to the antivortex in a pair to be along the x axis. The 
optimal coefficients correspond to Cy = (— 1) J , where 
the VA pairs orient themselves like dipoles in an anti- 
ferroelectric material. All topological dipoles are aligned 
along the x axis and anti-aligned along the y axis, as 
shown in Fig. |21 However, there is no easy axis as the 
system is rotationally invariant. The configuration where 
VA pairs are aligned as dipoles in a ferroelectric material 
is higher in energy. 

As the temperature is raised towards Tm, the impor- 
tance of vibrations and defects of the VA lattice increases 
and eventually causes its melting into a VA liquid state. 
In order to obtain the correct equations of motion, we 
have to recall that vortices and antivortices move perpen- 
dicularly to the applied force |2l| . Thus, their dynamics 
is Eulerian instead of Newtonian as in a lattice of atoms. 
This leads to the equation of motion 



fij-r- = — z x Fi, 
dt M 



(10) 



where F s ; = — V ri ninjG(ri — Tj) plays the role of 

the force. We can define displacement fields similar to 
the phonon problem in solids as u p (Ri) = r p (Rj) — R^, 
with the index p = V, A indicating a vortex (V) or an 
antivortex (A) sublattice, and R; labeling the equilib- 
rium lattice sites. This equation produces two modes 
in the long wavelength limit: a longitudinal (L) mode 
with frequency = c^q, with cj, = (ft/M)(7/,/&), 
and a tranverse (T) mode with frequency lot — exq, 
with ct — {h/M)(~fT /b), where ji, < jt are dimcn- 
sionless constants. The corresponding eigenvectors are 
T]l = uy + ua , and r\T = uy — . The harmonic action 
in diagonal form is 

Vi p § / 2 P 2 2 2 21 

Svh = ~tfT d r l a Li L ij + a L 2 L u + a Tl T i:j + a T2 T u \ , 

where < aA 4 < 1 are dimensionless elastic con- 
stants obtained from the interaction potential, and Ay = 
^[9ir]A,j + 9jT]A,i], with A = L,T. The VA lattice melting 
transition is then controlled by the Kosterlitz-Thouless- 
Halperin-Nelson- Young (KTHNY) mechanism [T3, E, 
[l9j. The melting temperature Tm is determined by the 
lowest energy mode (L) via the relation 



Tm = Ps 



«Li(«Li + 2a L2 ) 



(a Ll + a L2 ) 



(11) 
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FIG. 3: Square vortex-antivortex lattice in the strong cou- 
pling limit. 

together with Eqs. © and JSJ). The velocities of the VA 
lattice vibrational modes are related to the dimension- 
less elastic constants through = and 7t = 
Joer[j% The measurement of sound velocities can be use- 
ful to extract information about the short-range nature 
of the interaction forces through and . The values 
used in Fig. ^ correspond to Tm — 0-3p s {^, Ao, Tm)- 

In summary, we have studied the superfluid state of 
Fermi gases in a strongly confining one-dimensional op- 
tical lattice that suppresses quantum mechanical atom 
transfer (tunneling), making the Fermi system two- 
dimensional. In the s-wave channel, we have found that 
the Fermi gas undergoes a superfluid BKT transition and 
that vortex-antivortex (VA) bound states appear below 
3bkt during the evolution from weak to strong coupling. 
In the strong coupling limit, Tbkt is large and controlled 
by the Fermi energy of the atomic system, thus allowing 
for spontaneous VA pairs to be observed directly (with- 
out stirring the condensate). In addition, we have found 
that the ground state of the system evolves smoothly 
from weak to strong coupling, but supports a square VA 
lattice. Unlike the case of a lattice of identical vortices, 
which is triangular, a square VA lattice is favored, be- 
cause topological charge frustration pushes the triangu- 
lar VA configuration to higher energy. In addition, we 
have shown that the VA lattice melts via the KTHNY 
mechanism of dislocation mediated melting. Because the 
symmetry of the VA lattice is square, we expect that 
the melted lattice will not be hexatic, but quartic. The 
VA lattice might be detected in a trap release experi- 
ment without the need of rotating the Fermi superfluid. 
Experimental chances of this observation are higher in 
the strong coupling limit, where the melting tempera- 
ture Tm < ?bkt is also controlled by the Fermi energy. 



Although we have not yet studied in detail the possible 
quartic liquid VA state, there can be an additional tran- 
sition into an isotropic VA liquid before the BKT tran- 
sition is reached. In case experiments can be done using 
a velocity sensitive probe instead of a density probe as 
usual, not only the vortex cores would be seen, but also 
the direction of their rotation (vorticity) . Lastly, in addi- 
tion to the spin wave mode, we have found two VA lattice 
vibrational modes, one longitudinal and one transverse. 
All these modes disperse linearly, but the presence of a 
transverse mode is also a direct indication of the exis- 
tence of the VA lattice. 

We thank NSF (DMR 0304380) for financial support. 
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